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Abstract. We present X-ray spectral fits to a recently obtained Chandra grating spectrum of rj Carinae, one 
of the most massive and powerful stars in the Galaxy and which is strongly suspected to be a colliding wind 
binary system. Hydrodynamic models of colliding winds are used to generate synthetic X-ray spectra for a range 
of mass-loss rates and wind velocities. They are then fitted against newly acquired Chandra grating data. We find 
that due to the low velocity of the primary wind (s 500 kms' 1 ), most of the observed X-ray emission appears to 
arise from the shocked wind of the companion star. We use the duration of the lightcurve minimum to fix the wind 
momentum ratio at r, = 0.2. We are then able to obtain a good fit to the data by varying the mass-loss rate of the 
companion and the terminal velocity of its wind. We find that M 2 = 10 M© yr and v X2 ~ 3000 kms 
With observationally determined values of * 500-700 kms" 1 for the velocity of the primary wind, our fit implies 
a primary mass-loss rate of M , « 2.5 x 10" 4 M© yr' 1 . This value is smaller than commonly inferred, although 
we note that a lower mass-loss rate can reduce some of the problems noted by Hillier et al. (2001) when a va ue 
as high as 10~ 3 M© yr -1 is used. The wind parameters of the companion are indicative of a massive star which 
may or may not be evolved. The line strengths appear to show slightly sub-solar abundances, although this needs 
further confirmation. It also appears that the primary star is responsible for the nebula. 
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1. Introduction 

The superluminous star rj Carinae (HD 93308, HR 4210) 
continues to be extensively studied over a host of differ- 
ent wavelengths, yet remains intriguingly enigmatic. It is 
amongst the most unstable stars known. In the 1840s, and 
again in the 1890s, it underwent a series of giant outbursts 
(e.g. Viotti 1995) which ejected large masses of material 
into the surrounding medium. HST images of the result- 
ing bipolar nebula (e.g. Morse et al. 1998), known as the 
Homunculus, show it to be amongst the most spectacular 
in our Galaxy. The central object is now largely obscured 
by dust, and the cause of the outbursts and the nature of 
the underlying star (at outburst and today) remain spec- 
ulative. The source continues to show brightness fluctua- 
tions and emission-line variations. Further details of rj Car 
can be found in the review by Davidson & Humphreys 
(1997). 

In recent years evidence for binarity in this system 
has been accumulating. Damineli (1996) first noted a 
5.5 yr period in the variability of the He I 10830 A line. 
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Further photometric and radial velocity studies (Damineli 
et al. 1997, 2000), X-ray observations (Tsuboi et al. 1997; 
Corcoran et al. 2000 and references therein), and radio da- 
ta (Duncan et al. 1995, 1999) have supported the 5.5 yr 
period and the binary hypothesis (although, as noted 
later, recent HST observations do not match the 
ground-based radial velocity curve). A comparison of 
the abundances from the central object(s) and the com- 
position of the Homunculus nebula has also determined 
that there are at least two stars in this system (Lamers et 
al. 1998). 

rj Car is classified as a luminous blue variable (LBV). 
These are massive stars believed to be in a rapid and un- 
stable evolutionary phase in which many solar masses of 
material are ejected into the interstellar medium over a rel- 
atively short period of time (« 10 4 yr). LB Vs are regard- 
ed as a key phase in the evolution of massive stars, dur- 
ing which a transition into a Wolf-Rayet star occurs (e.g. 
Langer et al. 1994; Maeder &: Meynet 2001). Due to 
their rarity and complex nature however, we unfortunately 
still have no definitive theory for mass-loss during the LBV 
stage. The majority of the proposed mass-loss mecha- 
nisms are concerned with the importance of radiation pres- 
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sure within the outer envelope of the LBV, and for exam- 
ple utilize pulsational instabilities (e.g. Guzik et al. 1999), 
dynamical instabilities (e.g. Stothers k Chin 1993), or 
presuppose Eddington-like instabilities. The latter could 
arise from an enhancement in opacity as the star moves 
to lower temperatures (e.g. Lamers 1997), or from the in- 
fluence of rotation (e.g. Langer 1997; Zethson et al. 1999). 
Alternatively, the possibility that binarity plays a funda- 
mental role in explaining observed LBV outburst prop- 
erties has also been considered (Gallagher 1989), though 
most LBVs are not known binaries. Clearly, determining 
the wind and stellar properties of LBV stars is paramoun- 
t (see, for example, the discussions in Leitherer et 
al. 1994 and Nota et al. 1996). 

An important question is whether binarity influences 
the properties of LBVs (i.e. do LBVs in binaries evolve 
differently than single LBVs?). So while the presence of 
a companion can be exploited to help measure the mass 
of such stars, we must bear in mind that binary LBVs 
and single LBVs may be quite distinct objects (unless it 
is later discovered that all LBVs are actually binaries). 
Therefore, in order to use 77 Car to understand some of 
the defining LBV characteristics such as their extremely 
high mass-loss rates, we first need to determine beyond all 
doubt that 77 Car is in fact a binary, and then to determine 
the influence of the companion on the system. 

Investigations over the last few years have already 
helped to form a basic picture of 77 Car. The orbital 
parameters, although uncertain, indicate the presence of 
an early-type companion star, which will also have a 
powerful stellar wind. In such binaries, a region of hot 
shocked gas with temperatures in excess of 10 million K 
is created where the stellar winds collide (Prilutskii k 
Usov 1976; Cherepashchuk 1976). The wind-wind colli- 
sion (WWC) region is expected to contribute to the ob- 
served emission from this system, particularly at X-ray 
and radio wavelengths. Previous X-ray observations 
revealed extended soft emission from the nebula 
and strong, hard, highly absorbed, and variable e- 
mission closest to the star (Corcoran et al. 1995; 
Weis et al. 2001), in contrast to the emission char- 
acteristics from single stars, which are typical- 
ly softer, much less absorbed, substantially weak- 
er and relatively constant in intensity. Since 1996 
Feb 77 Car has been continuously monitored by RXTE in 
the 2-10 keV band (e.g. Corcoran et al. 2001a). The 
lightcurve (Fig. 1) contains remarkable detail showing a 
slow, almost linear, rise to maximum over a period of 
% 1 yr, followed by a rapid drop to approximately 1/6 
of the peak intensity for « 3 months, an almost as sharp 
rise to approximately 1/2 of the peak intensity level, and 
then almost constant intensity for « 3/4 of the proposed 
5.5 yr orbital period. The drop to minimum was successful- 
ly predicted from numerical models of the WWC (Pittard 
et al. 1998) before being actually observed. 

Small scale quasi-periodic outbursts in the X-ray 
lightcurve have also been detected (Corcoran et al. 1997). 
Estimates of changes in the timescale between succes- 
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sive flares as a function of phase were made by Davidson 
et al. (1998) for a variety of assumed orbital elements. 
RXTE X-ray observations obtained after the X-ray min- 
imum seem to show a lengthening of the flare timescale 
(Ishibashi et al. 1999), which indirectly support the bi- 
nary model. The latest published X-ray observation of 
77 Car is of a high resolution grating spectrum taken with 
the Chandra X-ray observatory (Corcoran et al. 2001b). 
Preliminary analysis has revealed the presence of strong 
forbidden line emission, which suggests that the density 
of the hot gas, n e < 10 14 cm- 3 . This can be contrasted 
with the newly published X-ray grating spectra of the sin- 
gle stars 0 1 Ori C (Schulz et al.2000), C Ori (Waldron k 
Cassinelli 2001), and £ Pup (Kahn et al. 2001), all of which 
have weak forbidden lines (indicative of either high den- 
sities or high UV flux near the line forming region). This 
lends further support to a wind-wind collision model. 

While it appears increasingly likely that 77 Car is in 
fact a binary system, using observations at lower energies 
to constrain the orbital elements and stellar parameters 
has proved difficult. Ground-based observations (for which 
good phase coverage exist) are hampered by poor spatial 
resolution and thus suffer contamination from strong neb- 
ular emission. High-spatial resolution spectra have been 
obtained by STIS but phase coverage is currently very 
limited. HST STIS observations at two different phases 
of the 5.52 year cycle (Davidson et al. 2000) did not con- 
firm the predicted variations in the radial velocity of the 
emission lines based on the ground-based radial velocity 
curve (Damineli et al.2000). If 77 Car is a binary, it is vi- 
tally important to determine the stellar parameters of the 
companion so that the effect of the companion on observa- 
tions can be understood and the correct stellar parameters 
of the primary can be derived. 


1.1. X-ray Emission from Colliding Winds 

The wealth of information contained in X-ray spectra of 
colliding wind binaries (e.g. the density, temperature, ve- 
locity, abundance, and distribution of the shocked gas in 
the wind collision region) has been a strong motivating 
force for observers and theorists alike in this field. Since 
the hot plasma in most colliding wind shocks is optically 
thin and collisionally ionized, and is generally assumed to 
be in collisional equilibrium, Raymond-Smith (Raymond 
k Smith 1977) or MEKAL (Mewe et al. 1995) spectral 
models are normally fitted to such data (e.g. Zhekov k 
Skinner 2000; Rauw et al.2000; Corcoran et al. 2001b). 
However, the multi-temperature, multi-density nature of 
the WWC region means that at best simple fits with 
one- or two-temperature Raymond-Smith models can on- 
ly characterize the broad properties of the emission. In 
this way one can estimate an ‘average’ temperature of the 
shocked gas, and an ‘average’ absorbing column, but lit- 
tle is learned of the underlying stellar wind parameters. 
At worst the application of one- or two-temperature mod- 
els to what is inherently multi-temperature emission can 
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Fie 1 Liehtcurve of r, Car observed with the RXTE satellite and phased to the 5.5 yr orbital period. A predicted hghtcurve 
(pfttard et al. 1998) from a numerical model of the wind-wind collision is overplotted. The two agree well, particularly the 
duration of the minimum. The rise from minimum is not in good agreement, but this is thought to be due to the limitations of 
modelling the wind collision in 2D. The rapid change in position angle of the stars through penastron passage leads to a greater 
absorbing column in fully 3D models (Pittard 2000). 


lead to spurious values of some of the fit parameters, e.g. 
abundances (cf. Strickland & Stevens 1998). 

Complex numerical hydrodynamical models have often 
been applied to gain insight into colliding wind system- 
s (e.g. Stevens et al. 1992, Pittard et al.1998). However, 
while undoubtedly useful, their interpretation can be d- 
ifficult, and to date there have been only two published 
papers where observed spectra are directly fitted with syn- 
thetic spectra from such models. In the pioneering work 
of Stevens et al. (1996), medium resolution ASCA spectra 
of the Wolf-Rayet binary 7 2 Velorum were fitted against a 
grid of synthetic spectra. In this fashion they were able to 
obtain direct estimates of the mass-loss rates and termi- 
nal velocities of the individual stellar winds. As mass-loss 
rates obtained from measures of radio flux or spectral line 
fits depend on a variety of untested assumptions, the im- 
portance of a new independent method to complement 
estimates from free- free radio or sub-mm observations, or 
from Ho: or UV spectral line fitting, cannot be stressed 
enough. Rates derived by Stevens et al. (1996) with this 
new method were significantly lower than the commonly 
accepted estimates for 7 2 Velorum based on radio obser- 
vations, but an indication of the future benefits of this 
method was realized when both sets of estimates were lat- 
er brought into agreement following a surprisingly large 
reduction in the distance to this star from Hipparcos da- 
ta 1 . We note that this method can also provide insights 
into the values of parameters which are otherwise virtu- 
ally impossible to estimate, such as the mass-loss rate of 


1 The thermal radio flux, S v oc M 4/3 D 2 (where D is the 

distance to the source) whereas the X-ray flux from an adia- 

batic wind collision is F v oc M 2 D~ 2 . Therefore M radio <x D 3 ^ 2 
whereas M xr ay oc D. If D is revised downwards, M radio de- 
creases faster than M xray • 


the companion, M 2 , or the characteristic ratio of the pre- 
shock electron and ion temperature (Zhekov & Skinner 
2000). The quality of recently available X-ray grating spec- 
tra now gives us access to important X-ray emission line 
diagnostics which should severely constrain models of the 
X-ray emission distribution. This means that stellar wind 
parameters can in principle be reliably estimated from 
analysis of X-ray grating spectra of colliding wind bina- 
ries. 

In addition to testing the binary hypothesis, the 
Chandra grating spectrum of r) Car (Corcoran et al. 2001b) 
provides the ideal opportunity to test the method devel- 
oped by Stevens et al. (1996) against a spectrum of much 
higher spectral resolution, and to pin down important 
physical parameters of the system. In this paper we there- 
fore fit the X-ray grating spectrum using a grid of colliding 
wind emission models to i) test the binary hypothesis, and 
ii) to attempt to obtain accurate estimates of the wind pa- 
rameters of each star. The fact that we are able to obtain 
good fits, with sensible model parameters, gives us further 
confidence in the binary hypothesis. We also find that un- 
like the UV and optical work where fits to the primary are 
made difficult by significant contamination from the com- 
panion star, instead the X-ray emission arises from the 
shocked wind of the companion and suffers essentially ze- 
ro contamination from the wind of the primary. Therefore 
the X-ray data uniquely samples parameters of the com- 
panion, in contrast to the optical analysis which probes 
the nature of the primary. In this sense our analysis is en- 
tirely complementary to the complex fits to the UV and 
optical HST spectrum of r/ Car by Hillier et al. (2001). Our 
analysis also provides us with a new estimate of the mass- 
loss rate of the primary star. For details of the Chandra 
observation and an initial analysis of the data the reader 
is referred to Corcoran et al. (2001b). In Section 2 we dis- 
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cuss the creation and variation of the model spectral grid; 
in Section 3 we describe the fit results; and in Section 4 
we summarize and conclude. 

2. The Synthetic Colliding Wind Spectra 

The method applied by Stevens et al. (1996) as discussed 
in Section 1 is as follows. We first calculated a whole host 
of hydrodynamical models with different stellar wind pa- 
rameters. We then generated a synthetic X-ray spectrum 
from each model. Finally, we fit the grid of synthetic spec- 
tra to the actual data. Best-fit values are then obtained for 
each parameter of interest (e.g. M i, ilf 2 , v OQl , v ^ etc. ). 
Specific points in relation to this method are highlighted 
below. 


greement with the 'expected value from our current un- 
derstanding of the orbit. We therefore fixed it at this val- 
ue for the rest of our analysis. However, we found that 
there were large uncertainties on the values of M x and 
v <»i; With the benefit of hind-sight this is not too sur- 
prising given the known slow speed of the primary wind 
(vod & 500 — 700 kms' 1 ), which means that the shocked 
primary wind is not a strong source of X-ray emission at 
hard energies. Therefore for our final grid we fixed the 
terminal velocity of the primary star at 500 kms' 1 and 
adjusted M i to obtain a desired value for the wind mo- 
mentum ratio, 7) = M 2 v 00 J{M 1 VoQl ), with M 2 and 
^002 as free parameters. 

The range of the free parameters used in our final grid 
is given in Table 1. To restrict the number of models to a 


2.1. Hydrodynamical Models of r\ Carinae 

The colliding wind models were calculated using VH-1 
(Blondin et al. 1990), a Lagrangian-remap version of the 
third-order accurate Piecewise Parabolic Method (PPM; 
Colella & Woodward 1984). The stellar winds are mod- 
elled as ideal gases with adiabatic index 7 = 5/3. 
Radiative cooling is included via the method of opera- 
tor splitting, and is calculated as an optically thin plasma 
in ionization equilibrium. The cooling curve for the tem- 
perature range 4.0 < log T < 9.0 was generated using 
the Raymond-Smith plasma code. Despite some evidence 
to the contrary (see Corcoran et al. 2001b), we have used 
solar abundances throughout this work, since it keeps the 
problem as simple as possible in this first investigation. 
We aim to model non-solar abundances in a later paper. 

The simulations were calculated assuming cylindri- 
cal symmetry of the wind collision zone - orbital effects 
are negligible at the phase of the Chandra observation 
(<A = 0.60). Grid sizes spanned a range from 300 x 300 to 
550 x 450 cells. Each grid cell was square and of constant 
size on an individual grid. The linear dimension of each 
cell was either 2.67 x 10 12 or 4.0 x 10 12 cm. The maximum 
distance from the axis of symmetry was 1.2 x 10 15 cm 
in all cases. At the phase of the observation, the orbital 
separation of the stars using the ephemeris of Corcoran et 
al. (2001a) is ss 4 x 10 14 cm. Thus there is either 100 or 150 
grid cells between the stars. The large separation means 
that the winds are likely to collide at very near their termi- 
nal velocities. We therefore assume that we can treat the 
winds as being instantaneously accelerated to their termi- 
nal velocities, and do not consider any radiative driving ef- 
fects. In comparison with the work on j 2 Velorum (Stevens 
et al. 1996), assuming terminal velocity winds and negli- 
gible binary rotation is more valid for 77 Car. We further 
assume that the winds are spherically symmetric. 

To compute a grid of synthetic spectra we initially 
varied four parameters (M u M 2 , Vo0i , Vo02 ) in the hy- 
drodynamic models. During our investigation we also in- 
cluded the separation of the stars, D, as an additional 
free parameter. From fits at this point it was found that 
D ^ 4 x 10 14 cm with a small uncertainty, in good a- 


manageable number the parameter steps are fairly coarse. 
In future papers we will use a much finer grid. The range 
in the value of 77 corresponds to either the wind of the 
primary (77 = 0.1) or the secondary (77 = 5) dominating, 
or the winds having equal momentum fluxes (77 = 1). The 
distance of the stagnation point from the centre of the 
primary star is given by 


r = 


1 + V^ 


D, 


( 1 ) 


and ranges from 1.2 x 10 14 cm to 3.0 x 10 14 cm. In the 
most extreme cases, this is still of order ten to a hundred 
times the radius of the stars, and justifies our use of termi- 
nal velocity winds. The kinetic power of the primary wind 
ranged from 4.7 x 10 34 ergs' 1 to 7.9 x 10 38 ergs' 1 , and 
from 7.1 x 10 35 ergs' 1 to 7.9 x 10 38 ergs' 1 for the sec- 
ondary wind. The combined kinetic power of the winds 
ranges from 7.6 x 10 35 ergs' 1 to 1.6 x 10 39 ergs' 1 . The 
half-opening angle of the contact discontinuity measured 
from the line between the secondary star and the shock 
apex ranges from « 50 - 120° (cf. Eichler & Usov 1993). 

The effect of radiative cooling can be quantified by 
the parameter x , the ratio of the cooling timescale to the 
dynamical timescale of the system. For shocked gas near 
the local minimum in the cooling curve at T rs 2 x 10 7 K, 


where i7 8 is the wind velocity in units of 1000 kms" 1 , 
d \2 is the distance to the contact discontinuity in unit- 
s of 10 12 cm, and M _ 7 is the mass-loss rate in unit- 
s of 10' 7 Mq yr' 1 (cf. Stevens et al. 1992). This equa- 
tion is valid for post-shock temperatures in the range 
10 7 < T < 10 82 K (680/ y/Ji <v< 2700/ ^j2 kms' 1 
where pm H is the average particle mass in grammes). For 
material with solar abundance (p = 0.6), this correspond- 
s to 0.9 < v s < 3.5. For winds with slower velocities, the 
shocked gas lies on the negative slope of the cooling curve. 
In this case the dependence of x on the velocity of the wind 
is steeper: 

^12 

a /_ 7 


(3) 
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Table 1. Our final grid of stellar wind parameters for the 
hydrodynamical models used to generate synthetic spectra. 
The models fix the terminal velocity of the primary wind at 
500 km s’*" 1 and the stellar separation at D ~ 4.0 x 10 14 cm. 
The mass-loss rate of the primary star is given by Af i = 
M 2 ^oo 2 /(500 rj). 


Parameter 

1 

2 

3 

n 

0.1 

1.0 

5.0 

M 2 (M 0 yr -1 ) 

nr 6 

10 -5 

10" 4 

V 002 (kms -1 ) 

1500 

3000 

5000 


In our grid of simulations, the shocked wind of the prima- 
ry star is almost always strongly radiative (xi.min 1>0) 
while the secondary’s is never strongly radiative (x 2 ,min = 
2.0). However, for some choices of the 3 parameters in 
Table 1, the primary’s wind can approach the point where 
it starts to become adiabatic (xi,max ~ 1*0). In Fig. 2 we 
show two examples of the hydrodynamic calculations used 
to calculate synthetic spectra. Due to the difference in ter- 
minal velocity of the two winds, the velocity shear at the 
wind interface always generates Kelvin-Helmholtz insta- 
bilities in our models. The rapid cooling of at least one 
of the two winds means that thin-shell instabilities also 
always occur in our models. The colliding winds region is 
most unstable when both shocked winds rapidly cool (cf. 
Stevens et al. 1992). 

The spectrum from each hydrodynamic model was av- 
eraged over 3 ‘snapshots’ each spaced by about 70 d. This 
is of order the wind dynamical timescale and serves to 
approximate a time-averaged spectrum. It is particularly 
important to adopt this approach for models where the 
wind collision region is very unstable since in these cases 
the flux can vary by greater than ±50% from one snapshot 
to the next. As already noted by Corcoran et al. (2001b) 
no variability is seen during the « 90 ksec exposure of the 
Chandra grating spectrum. This is to be expected, how- 
ever, since the dynamical timescale is much longer. Two 
examples of ‘time-averaged’ synthetic spectra are shown 
in Fig. 2, next to density plots of the corresponding hy- 
drodynamic calculation. 

The synthetic spectra were calculated over the 
energy range 1.26-10 keV (10-1.26 A respective- 
ly), and have a resolution of 0.005 A . The actual 
grating spectrum shows significant absorption be- 
low 1.5 keV and contains essentially no useful in- 
formation below our lower limit of 1.26 keV. The 
resolving power of our synthetic spectra is approx- 
imately twice as high as the grating data that we 
model. Use of the Raymond-Smith plasma code 
implicitly assumes that the plasma is in thermal e- 
quilibrium. This is generally true for colliding wind 
binary systems (see Pittard 2000). 
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2.2. Variations with rj, M 2 , and v CX)2 

In Fig. 3 and the following subsections we reveal how the 
shape and normalization of the calculated spectra depend 
on the various free parameters of the grid. 

2.2.1. Variation with t) 

Due to the low preshock velocity of the primary wind, the 
X-ray emission from the WWC is almost entirely from 
the shocked secondary wind. A greater fraction of the sec- 
ondary wind passes into the wind collision zone as the 
relative momentum of the primary wind is increased, so 
we expect the luminosity to increase as 77 decreases, as in- 
deed the top panel of Fig. 3 shows. An important point, 
however, is that the shape of the spectra do not seem de- 
pendent on the value of j], contrary to our expectations 
before starting this study. 

2.2.2. Variation with M 2 

As we increase the mass-loss rate of the secondary star in 
our models, the kinetic power and the density of its wind 
also increase. Therefore, for a given wind momentum ratio, 
the rate of conversion of kinetic energy to thermal energy 
increases, as does the density of the postshock gas. Thus 
there is more available energy to radiate and a greater ef- 
ficiency of radiation. Hence we should expect the X-ray 
luminosity to scale strongly with increasing M 2 - This is 
indeed what we see in the middle panel of Fig. 3, where an 
order of magnitude increase in Af 2 leads to an almost two 
order of magnitude increase in X-ray luminosity, indicat- 
ing that the calculations in this panel are in the adiabatic 
regime (L x oc M 2 ; see Stevens et al. 1992). However, there 
is no evidence for a softening of the spectrum with increas- 
ing M 2 - This is again contrary to our initial expectations 
(note, however, the discussion in the next section). 

2.2.3. Variation with Voo 2 

If the preshock velocity of a wind is increased, the post- 
shock temperature increases (T oc v 2 ). Therefore we ex- 
pect the model spectra to harden as the preshock velocity 
of the secondary wind is increased. Because the kinetic 
power of the wind increases we also expect the luminosi- 
ty to increase. The variation of the spectrum with v ^ is 
shown in the bottom panel of Fig. 3. Here we find that 
our expectations are only partially met. The overall lu- 
minosity does increase with Uoo 2 > and the spectrum does 
harden between Voo 2 = 1500-3000 kms -1 , but it softens 
between v ^ = 3000-5000 krns -1 . An investigation into 
this surprising finding offers the following explanation. 

For fixed rj and M 2 , and variable v^, the density of 
the wind varies as p oc M /v^. If v. ^ is increased p de- 
creases, which reduces the radiative efficiency of the hot 
shocked gas near the stagnation point ( E oc n 2 ), which 
acts to slightly suppress the hard flux. A plot of the aver- 
age density against T for v X2 = 1500,3000,5000 kms 
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r*±? y t° A r m>C “T a " d theoretical s P ectra ° f the colliding winds in r, Car. The top panels show the results from 
w ere the momentum of each wind is equal (17 = 1 ), whereas in the bottom panels the primary wind has a momentum 10 

rthe^ft a a!h T ™ * = ° ^ *“ b0th * density plot (g Cm ^ from the hydrodynamical grid is shown 

on the left and the resulting intrinsic X-ray spectrum shown on the right (no absorption). On each hydrodynamical grid the 

primary star ,s located at (0,0) and the secondary star is at (4.0 x 10 14 ,0). The shape and structure of the shocked region varies 

' th m Tf nl balance between the wmds and the respective value of * for each shocked wind. The full parameters used 

in these models were- imnor nanolc K.f . in — 3 * * — 1 « ^ — 4 _t , " useu 


in these models were: upper panels - A/i = 10 

lower panels - M i = 3 x 10 4 M© yr _1 , M 2 = iu ~ m© yr * , v 00l — „„ 
shape of the X-ray spectrum are a direct consequence of these parameters. 


3 Moyr" 1 , M 2 = 10~ 4 M q yr -1 , 
10 5 Moyr -1 ,!;^, = 500 kins -1 


v <x t — 500 kms 1 , Voo 2 - 5000 kms" 
v oo 2 ~ 1500 kms 1 . The luminosity and 


also reveals that for 10 6 3 < T < 10 7 4 the average densi- 
ty is higher when Voo 3 = 5000 kms -1 than when v 0O2 — 
3000 kms 1 , which tends to keep up the soft emission 
relative to the hard emission. This appears to be because 
there is some mixing at the interface between the prima- 
ry and secondary wind, and the shocked primary wind is 
denser when v OQ2 = 5000 kms -1 {M \ = 3 x 10 — 4 , 6 x 
10- 4 , lO- 3 Moyr' 1 for v X2 = 1500,3000,5000 kms’- 1 ). 
The mixed gas tends to be at these intermediate temper- 
atures. 

The overall net effect is that the emission from the hot 
gas is suppressed relative to the cooler gas, which softens 
the spectrum from the model with v ^ = 5000 kms -1 
relative to the model with v X2 = 3000 kras -1 . These in- 
ferences are supported by spectra calculated with = 
5000 kms -1 and tj = 0.1, and variable M 2 (Fig. 4). This 
shows that the spectrum hardens as A/ 2 increases from 
10 6 M w yr _1 to 10 -4 M 0 yr -1 (the spectral index 


a, where F„ oc 1/", increases by about 0.1). This is 
what we expect from our prior inferences since a higher 
M 2 helps to make postshock gas with a preshock velocity 
°f ^oo 2 = 5000 kms -1 somewhat more efficient at radi- 
ating. In summary, we find that varying M 2 can actually 
have a small effect on the shape of the spectrum in some 
parts of parameter space, as well as the more obvious and 
much greater effect on the luminosity. The power of M 2 
on affecting the spectral shape is, however, much less than 
that of Uoo 2 , which is the primary influence. 

2.2.4. Summary of Spectral Dependence on rj, M 2 , 

and Wqoj 

For the models investigated in this paper, we find that 
the spectral shape is i) insensitive to the value of 77 , ii) 
can harden slightly as M 2 increases (in other parts of pa- 
rameter space, e.g. strongly radiative shocks, it may soften 





Fie. 3. Variation in the spectral shape and normalization of the theoretical colliding winds spectra in 77 Cannae. The top panel 
shows the variation with 77 (M 2 = 10 ' 6 Moyr ' 1 , t >« 2 = 3000 km s' 1 ). The middle panel shows the variation with M 2 
(t) = 5.0, Doo 2 = 1500 kms -1 ). The bottom panel shows the variation with 11002 (77 = 0.1, M 2 - 10 M© yr ). 


with increasing M 2 ), and iii) generally hardens with in- 
creasing U 002 , but can sometimes soften. The luminosity 
increases with i) decreasing 7 ? (since 77 Car is unusual in 
the sense that there is no discernible contribution to the 
X-rays from the shocked primary wind, bar mixing), ii) 
increasing M 2 , and iii) increasing i>oo 2 - 


2.3. The Spectral Grid 

The synthetic spectra were input into a FITS-style ‘Table’ 
model for use with the XSPEC data-analysis package, as 
specified by Arnaud (1999). The fitting procedure with- 
in XSPEC is able to interpolate between the parameter 
values on the spectral grid so that best-fit values of M 2 
etc. are obtained from fitting the Chandra data. In the 
2-10 keV band, the intrinsic luminosity on our grid 
ranges from 1.8 x 10 31 ergs' 1 to 1.6 x 10 37 ergs' 1 , 
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Time (years) 

Fig. 5. The duration of the X-ray minimum predicted from 
simple models. The wind momentum ratio in the three models 
shown is ij = 0.2, TO, 5.0 (solid, dots, dashes), which gives a 
duration of 92,33, 16 d respectively. The observed duration is 
~ 100 d, so the data is best approximated bv r/ = 0.2. 


4- Variation in the spectral shape and normalization for 
a very fast secondary wind (voo 2 = 5000 kms -1 ) and a fixed 
wind momentum ratio (t) = 0.1). As M 2 increases the spec- 
trum becomes slightly harder. 

which brackets the observationally determined value of 
4 x 10 34 ergs -1 (Corcoran et al. 2001b). 

3. Fits to the Chandra grating spectrum 

We compared the model intrinsic colliding wind spectra 
with the observed grating spectrum allowing the mass- 
loss rate and wind terminal velocity of the companion to 
vary. We in addition included absorption from a cool over- 
lying medium to simulate the combined absorption of the 
circumstellar (nebular) material in the Homunculus, the 
ISM, and the wind from the companion (all of which may 
contribute to the X-ray absorption at some level). Since 
the stellar wind absorption depends on the companion’s 
mass-loss rate (at the phase of the Chandra observation) 
for consistency we should model this component of the 
absorption directly, but this is beyond the scope of the 
present paper. Because of the low number of counts per 
spectral bin in the high-resolution grating spectrum, and 
because background does not contribute significantly to 
the observed spectrum at energies above 2 keV, we fit the 
gross spectrum using the C-statistic (Cash 1979), which 
is appropriate for Poisson-distributed data. 

The fact that the shapes of our model spectra are not 
dependent on the value of t) introduces a degeneracy into 
our grid as both Tj and M 2 primarily influence the nor- 
malization of the spectra. This means that various com- 
binations of T) and M 2 can provide similar quality fit- 
s to the data. However, it is possible to obtain an esti- 
mate for Tj from the duration of the lightcurve minimum 
(~ 100 d). To do this we have constructed a simple model 
of the observed X-ray emission. We assume that the in- 
trinsic flux varies as l/£>, and that the absorption from 


the shock apex along the line of sight is negligible when 
viewing through the less dense secondary wind, but total 
when viewing through the very dense primary wind. We 
use the latest orbital parameters (Corcoran et al. 2001a) 
and assume that the skew angle of the shock cone from 
the orbital velocity, w or b, is 6 ~ arctan(u orb / Uoo J. Fig. 5 
shows the results, scaled so that L x = 1 at periastron. The 
skew of the shock breaks the symmetry of the observed e- 
mission so that the post-minimum flux is lower than the 
pre-minimum flux, as observed. The duration of the min- 
imum decreases with increasing ij (for 7/ = 0.1, 0.2, 0.3 the 
duration is 133, 92, 71 d), and is best matched by rj » 0.2 2 . 
Therefore we adopt this value for the rest of our analysis. 

We attempted to fit the grid of synthetic colliding 
wind spectra to the observed MEG +1 order spectrum. 
We also constrained the spectrum normalization to lie be- 
tween 0.8 and 1.2, which fixes the distance to the star 
between 2300 and 1900 pc, close to the canonical distance 
of 2100 pc. In Fig. 6 we show the best fit interpolated 
spectrum from our grid in the range 1.5 - 7.5 keV. Fig. 7 
shows a close-up of the fit to the Si and S emission lines be- 
tween 1.8-2.8 keV. Both figures demonstrate the excellent 
agreement found between the data and the models and 
highlight the progress made in both the observational and 
theoretical study of colliding winds over the last 10 years. 
The continuum shape and level from the models matches 
closely that of the data. All strong lines (and most weak 
ones) which appear in the observed spectrum are matched 

Incidently, the lightcurve minimum computed with the 2D 
hydrodynamical model in Pittard et al. 1998 (redisplayed in 
Fig. 1 at the beginning of this paper) also matches the du- 
ration of minimum very well (although not the post-minimum 
flux because the skew of the shock could not be computed with 
a 2D calculation). Although different mass-loss rates and a s- 
lightly less eccentric orbit were used, the wind momentum ratio 
adopted was also Tj = 0.2, which provides a certain robustness 
for this value. 
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Table 2. The best fit parameters from the fit of the spec- 
tral grid to the grating data. We fix the terminal velocity of 
the primary wind at 500 km s 1 , the stellar separation at 
£) = 4.0 x 10 14 cm, and the wind momentum ratio (ij) at 
0.2, and allow M 2 , Voo 2 , Nh and the normalization to vary. 
Solar abundances were assumed, and a cold external absorb- 
ing column was applied. The resultant normalization is very 
close to 1.0, which indicates that the model and the values of 
D and rj are sensible, and also that the assumed distance to 
r) Cm is approximately correct. The implied mass-loss rate of 
the primary star from this fit is M i = 2.5 x 10 M© yr 

Parameter Value 


V 

0.2 (fixed) 

M 2 ( M© yr -1 ) 

0.98^08 x 10-5 

1^002 ( S ) 

3000^340 

Nh (cm- 2 ) 

7.68±o;ir x 10 22 

Normalization 

1.16 


also in the best fit. This means that the temperature dis- 
tribution seen in the grating spectrum is matched by the 
model. However most of the strong lines are significant- 
ly overpredicted by the model fit. This discrepancy is in 
the opposite sense to the fit results from simple one- and 
two-temperature Raymond-Smith models, and is perhaps 
revealing that our assumpt ion of solar abundances needs 
to be modified. The best fit parameters are summarized 
in Table 2. 

4. Conclusions 

We have presented the first direct fits to the X-ray spec- 
trum of r) Car using a grid of spectra calculated from hy- 
drodynamical models of the wind-wind collision. This is 
only the second time that this technique has been applied 
to a colliding winds X-ray spectrum and the first time 
with a high quality grating spectrum. 

We find that the colliding wind emission model nat- 
urally provides for the range of ionization seen in the e- 
mission line grating spectrum for reasonable values of the 
wind parameters. Due to the low velocity of the primary 
wind (« 500kms _1 ), most of the observed X-ray emission 
arises from the shocked wind of the companion star. We 
find it difficult therefore to fit both M i and v OQl as free 
parameters. However, the duration of the observed X-ray 
minimum can be used to estimate the wind momentum 
ratio of the stars, q. With t) fixed at 0.2, and M 2> an< i 
Uoo 2 as free parameters, we are able to obtain a good fit 
to the data. 

We find that the mass-loss rate of the companion is 
M 2 ~ 10 -5 M© yr -1 and the terminal velocity of it- 
s wind is Uoo 2 « 3000 km s' 1 . These values suggest 
that the companion is probably of type 06 or earlier 
(cf. Howarth & Prinja 1989), or is possibly a WR star. 
The velocity of the primary wind has been determined 
to lie in the range 500 - 700 kins -1 (e.g. Hillier et 
al.2001). Hence our fit implies a primary mass-loss rate 
of M i « 2.5 x 10 -4 M© yr -1 . From the uncertainty in 


the value of r/, and the interpolation on our grid, we es- 
timate the uncertainty on our derived value for M i as 
approximately a factor of 2. 

Our best-fit estimate of M i is smaller than typical- 
ly inferred (cf. Davidson & Humphreys 1997; Hillier et 
al.2001) However, we note that a lower mass-loss rate 
for the primary star can reduce some of the problems 
noted by Hillier et al. (2001) who fitted a value as high 
as 10 -3 M©yr -1 . In particular, the models of Hillier 
et al. (2001) suffered from absorption components that 
were too strong and electron-scattering wings which were 
overestimated. Both indicate that their chosen mass-loss 
rate is too high. Paradoxically both the Ha and H/3 e- 
mission lines were weaker than observed, indicating that 
their mass-loss rate is too low. However, it is well known 
that the wind collision zone can be a strong source of 
emission lines (e.g. HD 5980, Moffat et al. 1998), which 
would resolve this problem. Their inferred minimum col- 
umn density is also larger than the observed X-ray value, 
again implying an overestimate of Mi- To increase the 
mass-loss rate of the primary towards the value estimat- 
ed by Hillier et al. (2001), we require either a reduction 
in the wind velocity to ~ 100 km s 4 , which in turn is 
in conflict with previous estimates, or a reduction in the 
wind momentum ratio to r) < 0.1, which is in conflict with 
the observed duration of the X-ray minimum. Finally, it 
is worth noting that our results indicate a value for M i 
which is closer to the value inferred for the Pistol star 
(M « 4 x 10 -5 Mg yr -1 ; Figer et al. 1998), an extreme 
early-type star with similarities to t] Car. 

Since current observations and theoretical modelling 
of the optical spectrum are unable to determine the ef- 
fective temperature and the stellar radius of the primary 
without first determining M \ (cf. Hillier et al.2001), our 
independent estimate may prove to be extremely useful 
in this regard. It will be interesting to see if the results 
from this paper are consistent with future X-ray obser- 
vations, and whether estimates of M x from observations 
at X-ray and other wavelengths can be reconciled. Future 
X-ray grating observations should also help us to fix the 
value of the wind momentum ratio more accurately. 

As the secondary wind dominates the X-ray spec- 
trum, and its terminal velocity appears to be high (« 
3000 km s _1 ), we should expect to see signs of Doppler 
broadening and shifts in the line profiles. While there is 
little evidence for this in the current spectrum, other or- 
bital phases may be more favourable in this regard. The 
addition of Doppler effects has already been incorporated 
in modelled X-ray spectra (Pittard et al.in preparation), 
and should provide further information on wind velocities 
and the structure of the wind-wind collision region. 

The over-prediction of the X-ray lines in our models 
indicates that the companion has sub-solar abundances, 
which favours an O-type over a VVR classification, al- 
though we would need to perform a more detailed analysis 
to confirm this possibility. As the primary has slightly en- 
hanced abundances of C and N compared to solar, this 
perhaps suggests that to date there has been no mass ex- 
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data and folded model 



channel energy (keV) 


r= I'9 B ? 7 he s P ectrum - The fit parameters were M 2 = 0.98 x 1<T 5 M Q yr~ l , v X2 = 3000 kras' 1 

If, = 2 ' 5 x H ln -/'L X vr -i T ’ an f ? °T a11 normallzation wa s 1.16. If we take Wool = 600 kms' 1 , this implies that 
o y r -As seen from the figure, the continuum level is very well fitted by the interpolated model. 


change between the stars. Lamers et al.(1998) suggested 
that the star which dominates the UV G HRS spectrum 
is not the star which ejected the nebula since the abun- 
dances in the GHRS spectrum are not as evolved as the 
abundances in the nebula (which are indicative of CNO- 
cycle products). As the UV bright source is probably the 
companion (Hillier et al. 2001) this indicates that it was 
the primary which ejected the nebula. There are also some 
caveats about the analysis in Lamers et al. (1998) since 
strong C lines (which Lamers et al.took to indicate nor- 
mal CNO abundances in the GHRS spectrum) also appear 
in stars known to be deficient in C (see the discussion in 
Hillier et al. 2001). It is also interesting to note that the 
UV spectrum brightened in 1999.1 vs. 1998.2, which is 
suggestive of an eclipse of the UV source near periastron 
in 1998. In conclusion, the high energy photons (UV and 
X-ray) seem to be telling us about the companion. 

We emphasize that contrary to the vast majority of col- 
liding wind systems, our X-ray analysis of tj Car primarily 
probes the conditions of the shocked wind of the compan- 
ion. X-ray observations of i / Car are therefore unique in 
this regard since at other wavelengths (with the possible 
exception of the far UV) the wind of the primary dom- 


inates the observed phenomena. While our analysis has 
for the first time provided a direct estimate of the wind 
parameters of the companion star, relating these to the 
stellar parameters (mass, radius, luminosity) of the com- 
panion star requires more work. It is anticipated that the 
continued multi-wavelength study of 77 Car through and 
beyond the next periastron passage will further reveal the 
hidden secrets of this most enigmatic system. 
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